AEROSPACE  REPORT  NO. 
TR-2006(8565)-3 


THE  AEROSPACE 

CORPORATION 

egundo,  California 


Measurements  and  Modeling  of  SiCl4 
Combustion  in  a  Low-Pressure  H2/02  Flame 


10  November  2006 


Prepared  by 

T.  MOORE,  B.  BRADY,  and  L.  R.  MARTIN 
Space  Materials  Laboratory 
Laboratory  Operations 


Prepared  for 

SPACE  AND  MISSILE  SYSTEMS  CENTER 
AIR  FORCE  SPACE  COMMAND 
483  N.  Aviation  Blvd. 

El  Segundo,  CA  90245-2808 


20070423141 

Engineering  and  Technology  Group 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


This  report  was  submitted  by  The  Aerospace  Corporation,  El  Segundo,  CA  90245-4691,  under  Con¬ 
tract  No.  FA8802-04-C-0001  with  the  Space  and  Missile  Systems  Center,  483  N.  Aviation  Blvd.,  El 
Segundo,  CA  90245.  It  was  reviewed  and  approved  for  The  Aerospace  Corporation  by  G.  F. 
Hawkins,  Principal  Director,  Space  Materials  Laboratory;  and  D.  C.  Marvin,  Principal  Director, 
Office  of  Research  &  Technology  Applications.  Michael  Zambrana  was  the  project  officer  for  the 
Mission-Oriented  Investigation  and  Experimentation  (MOIE)  program. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  (PAS)  and  is  releasable  to  the  National 
Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be  available  to  the  general  public,  including 
foreign  nationals. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication.  Publication  of  this  report 
does  not  constitute  Air  Force  approval  of  the  report's  findings  or  conclusions.  It  is  published  only 
for  the  exchange  and  stimulation  of  ideas. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


ubltc  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  this  coflection  of  information.  Sena  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  information  Operations  and 
Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person 
shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM 


TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  ( DD-MM-YYYY) 

10-09-2006 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  TYPE 


Measurements  and  Modeling  of  SiCl4  Combustion 
in  a  Low-Pressure  H2/02  Flame 


3.  DATES  COVERED  (From  -  To) 


5a.  CONTRACT  NUMBER 

FA8802-04-C-000 1 


Sb.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 


5d.  PROJECT  NUMBER 


T.  Moore,  B.  Brady,  and  L.  R.  Martin 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

The  Aerospace  Corporation 
Laboratory  Operations 
El  Segundo,  CA  90245-4691 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Space  and  Missile  Systems  Center 
Air  Force  Space  Command 
483  N.  Aviation  Blvd. 

El  Segundo,  CA  90245 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


TR-2006(8565)-3 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

SMC 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

Laser-Induced  Fluorescence  (LIF)  was  used  to  measure  temperature  and  OH  concentration  profiles  as  a  function  of 
distance  from  a  McKenna-style  burner  in  premixed,  one-dimensional,  low-pressure  H2/02/Ar  SiCl4-doped  flames.  The 
addition  of  SiCl4  was  shown  to  affect  the  flame  temperature  and  OH  concentration  profiles.  A  gas-phase  chemical 
kinetics  mechanism  for  the  combustion  of  SiCLj  in  an  H2/02/Ar  flame  was  proposed,  and  experimental  results  were 
compared  with  predictions  for  a  premixed,  one-dimensional  laminar  flame  model  based  on  CHEMKIN.  The  low- 
pressure  flame  data  are  sensitive  to  the  overall  kinetics  of  the  mechanism.  In  order  to  obtain  the  best  fit  to  the 
observed  data  for  all  flame  configurations,  we  had  to  modify  six  different  rates  from  our  original  estimates.  None  of 
the  modified  rates  are  well  known  for  the  temperature  regime  of  our  flame.  Particle  formation  and  surface  chemistry 
were  not  taken  into  account. 


15.  SUBJECT  TERMS 

Low-pressure  flame,  SiC14,  Silicon  tetrachloride,  Combustion  synthesis.  Flame  synthesis.  Kinetics  mechanism 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 


b.  ABSTRACT 


UNCLASSIFIED  UNCLASSIFIED 


C.  THIS  PAGE 


UNCLASSIFIED 


17.  LIMITATION 

18.  NUMBER 

OF  ABSTRACT 

OF  PAGES 

24 

RESPONSIBLE  PERSON 

Teresa  Moore 


19b.  TELEPHONE  NUMBER 
(include  area  code) 
(310)336-2486 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


Contents 


1 .  Introduction .  1 

2.  Experiment .  3 

3.  Mechanism  and  Model .  7 

4.  Results  and  Analysis .  13 

5.  Discussion .  17 

5.1  Reaction  Rates .  17 

5.2  Temperature .  19 

5.3  OH  concentration .  19 

6.  Conclusion . .  21 

References .  23 

Figures 

1 .  Schematic  of  low-pressure  flame  apparatus .  3 

2.  Fluorescence  data  collected  1 .2  cm  from  the  burner  for  five  data  runs .  13 

3.  Temperature  profiles  for  hydrogen/oxygen  flames  with  and  without  SiCl4 .  14 

4.  OH  concentration  profiles  for  hydrogen/oxygen  flames  with  and  without  SiCl4 .  15 


Tables 


1.  Silicon  Tetrachloride  Proposed  Combustion  Mechanism  (k  =  ATbe_E/RT) .  7 

2.  Transport  Database  for  Species  from  SiCl4  Combustion.  ()  =  estimated  quantity .  9 

3.  Thermodynamic  Data  for  the  Silicon  Species  in  Mechanism .  10 

4.  CHEMKIN  Sensitivities  of  OH  Concentration  to  the  Reaction  Rates 

at  Different  Distances  from  the  Burner .  17 

5.  Reactions  modified  in  order  to  improve  the  model  fit  to  our  data .  18 


IV 


1.  Introduction 


Formation  of  high-purity  powders  and  glasses  via  flame  aerosol  processes  is  important  for  many 
industries,  including  fiber  optics,  ceramics,  and  semiconductors.  Although  many  studies1,2  have  been 
carried  out  on  the  formation  of  silicon  dioxide  (Si02)by  burning  various  silicon-containing  precur¬ 
sors,  the  fundamental  chemical  kinetics  leading  to  particle  formation  are  not  well  understood.  Mod¬ 
eling  the  gas-phase  chemical  kinetics  is  difficult  due  to  the  paucity  of  fundamental  data,1  including 
the  primary  reaction  pathways  and  the  intermediate  species.  Silane  combustion  has  been  the  most 
studied;  several  detailed  chemical  kinetics  mechanisms  for  silane  combustion  varying  in  complexity 
and  physical  processes  considered  have  been  proposed. 3,4,5  Combustion  of  other  silicon  precursors  is 
less  well  understood. 

Many  of  the  experimental  and  modeling  efforts  on  the  formation  of  Si02  in  a  flame  have  focused  on 
determining  final  particle  size  and  morphology.  For  recent  reviews  on  the  formation  of  oxides  in 
flames,  see  Wooldridge1  and  Pratsinis.2  The  physical  flame  configuration  and  temperature  profile6 
determine  the  extent  to  which  precursor  identity  and  kinetics  affect  the  final  particle  size,  morphol¬ 
ogy,  and  composition.  Several  studies  have  shown  that  addition  of  a  silicon-containing  dopant  to  dif¬ 
fusion  and  premixed  flames  can  significantly  affect  flame  characteristics,  including  temperature  and 
species  concentrations.7,8,9,10  Because  temperature6  is  an  important  influence  on  the  final  particle 
characteristics  and  the  kinetics  affect  the  flame  temperature,  understanding  chemical  kinetics  of  the 
flame  is  important  for  understanding  and  predicting  final  particle  characteristics. 

There  are  several  models  that  deal  with  particle  formation,  but  to  date,  no  detailed  chemical  kinetics 
mechanisms  for  the  reaction  of  silicon  tetrachloride  (SiCU)  in  hydrogen/oxygen  flames  have  been 
published.  A  simple  chemical  kinetics  mechanism  based  on  an  initial  reaction  of  SiCU  +  H20  was 
proposed  by  Hannebauer  and  Menzel.1 1  SiCU  combustion  is  often  modeled  by  simplifying  the  kinet¬ 
ics  to  the  overall  oxidation  (SiCl4  +  02  -»  Si02  +  2  Cl2)  or  hydrolysis  (SiCU  +  H20  -»  Si02  +  4HC1) 
reactions,2  or  by  using  several  rate-determining  steps12  to  model  the  flame  chemistry. 

Various  less  complex  systems  involving  SiCU  +  02  or  SiCU  +  H2  reactions  for  lower-temperature 
chemical  vapor  deposition  applications  have  been  modeled,13,14,15  but  for  the  most  part,  these  models 
either  include  a  very  simplified  kinetics  scheme16,17  or  dismiss  kinetics  all  together.18  One 
model  ’  ■  proposed  a  low-temperature  SiCU  +  02  mechanism  in  which  large  chlorosiloxane 
chains  and  cyclic  compounds  are  formed,  progressively  becoming  more  and  more  like  Si02  with 
time.  However,  none  of  these  models  represents  the  conditions  in  a  flame  environment. 

While  these  approaches  to  modeling  can  yield  valuable  information,  a  more  detailed  model  of  the 
flame  chemical  kinetics  is  necessary  to  fully  understand  the  initial  particle  formation  process  and  how 
it  relates  to  the  final  particle  properties.  We  have  proposed  a  gas-phase  chemical  kinetics  mechanism 
for  the  combustion  of  SiCl4  in  an  H2/02/Ar  flame.  In  this  report,  we  report  measurements  of  tem¬ 
perature  and  OH  concentration  profiles  for  a  premixed,  one-dimensional,  low-pressure  flat  H2/02/Ar 
flame  doped  with  SiCl4.  We  compare  these  data  to  modeling  results  using  our  mechanism  and  the 
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CHEMKIN  modeling  software. t231,[241  The  goal  of  this  experiment  was  to  test  the  overall  reaction 
kinetics  of  our  model  in  a  one-dimensional,  premixed  flat  flame  where  the  role  of  fluid  mechanics 
was  simplified. 

Heterogeneous  chemistry  and  particle  formation  were  not  included  in  the  model.  This  is  not  because 
these  processes  are  negligible  or  unimportant,  but  because  there  is  no  good  treatment  for  modeling 

") S  7ft  97 

their  impact  on  a  low-pressure  flame.  There  are  good  modeling  treatments  for  these  processes,  ’  ’ 
including  the  CHEMKIN  package,  but  models  that  include  heterogeneous  processes  do  not  include 
the  transport  mechanisms  that  are  important  in  this  flame.  We  chose  the  current  model  because  it  is 
the  standard  tool  for  modeling  low-pressure  flames. 

We  used  a  low-pressure  flame  technique28  because  the  homogeneous  flame  reaction  zone  and  the 
particle  formation  region  is  significantly  elongated  compared  to  atmospheric-pressure  flames.  This 
increases  the  spatial  resolution  with  which  we  can  probe  the  temperature  and  species  concentrations. 
Our  main  diagnostic  technique  for  measuring  radical  species  concentrations  was  Laser-Induced 
Fluorescence  (LIF).29 
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2.  Experiment 


Figure  1  shows  a  diagram  of  our  low-pressure  flame  apparatus.  A  low-pressure  porous  plug  flat- 
flame  burner  was  mounted  on  a  three-axis  translator  in  a  vertical  flow  configuration.  The  chamber 
pressure  was  maintained  at  30+0.5  Torr;  it  was  servo-controlled  by  using  a  100  Torr  capacitance 
manometer  to  provide  feedback  to  a  throttle  valve  on  the  vacuum  pump  line.  The  mechanical  pump 
on  the  chamber  was  purged  with  nitrogen  and  equipped  with  mechanical  traps  to  capture  particles. 

The  combustion  gases  were  mixed  in  a  small  manifold  before  entering  the  burner.  The  gas  flows 
were  regulated  by  computer-controlled  mass  flow  controllers.  The  burner  outer  shield  flow  was  run 
at  2.7  standard  liters  per  minute  (SLM)  of  argon  to  help  keep  the  premixed  combustion  zone  one¬ 
dimensional.  The  total  flow  rate  through  the  burner  was  kept  constant  at  2.25  SLM. 


Three  basic  flame  configurations  were  studied:  an  undiluted  stoichiometric  H2/O2  flame,  a  diluted 
stoichiometric  H2/O2  flame  with  35%  of  the  total  flow  as  argon  diluent,  and  a  fuel-rich  undiluted 
flame  with  O  =  1 .2.  These  flames  were  studied  with  and  without  SiCl4  precursor  at  concentrations  of 
0.5%  of  the  H2+  O2  flow.  For  the  experiments  doped  with  SiCl4,  the  oxygen  flow  exiting  the  flow 
meter  was  bubbled  through  a  Teflon-coated  stainless-steel  cell  containing  SiCl4  before  being  mixed 
with  the  other  gases.  The  cell  was  temperature  controlled  to  within  ±0.005  °C.  SiCl4  concentration 
was  determined  from  the  vapor  pressure  and  total  pressure  measurements,  which  were  carried  out 
with  capacitance  manometers.  The  SiCl4  cell  was  weighed  before  and  after  each  experiment  and 
compared  to  the  calculated  mass  loss  to  confirm  the  calculated  SiCl4flow  rates. 
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Figure  1 .  Schematic  of  low-pressure  flame  apparatus. 


3 


Temperature  and  OH  concentration  profiles  were  obtained  using  a  standard  laser-induced  fluores¬ 
cence  (L1F)  technique  modified  for  use  in  a  particle-forming  combustion  environment.9  The  laser 
system  consisted  of  a  dye  laser  pumped  by  a  Nd:YAG  laser  operating  at  532  nm  and  30  Hz.  The 
resulting  output  was  doubled  with  a  KDP  crystal  and  then  irised  to  ~0.6  mm  diameter  before  being 
sent  into  the  vacuum  chamber.  The  fluorescence  was  collected  perpendicular  to  the  laser  axis  with  a 
lens  filtered  by  308±1 5-nm  bandpass  filters  and  then  focused  onto  a  photomultiplier  tube.  Signals 
were  processed  by  a  boxcar  integrator  triggered  by  the  Nd:YAG  laser.  A  short  (15  ns)  gate  coinci¬ 
dent  with  the  laser  pulse  was  used  in  order  to  reduce  quenching  effects  on  the  signal.  Power  meters 
recorded  the  power  both  before  and  after  the  chamber,  allowing  the  fluorescence  signal  to  be  cor¬ 
rected  for  fluctuations  in  laser  intensity.  Typical  laser  energy  used  was  0.3-2.0  pJ/pulse.  Experi¬ 
ments  were  performed  to  verify  that  these  energies  were  low  enough  to  avoid  saturation  of  the  OH 
transition.  The  optics  were  purged  with  argon  to  avoid  particle  deposition. 

For  temperature  and  OH  concentration  measurements,  select  transitions  in  the  (1  «-0)  band  of  the 
A  *-X system  were  excited  and  the  broadband  fluorescence  around  3 1 0  nm  was  observed.  For  temp¬ 
erature  measurements,  the  relative  intensities  of  the  P|(2)  and  Ri(14)  lines  near  282.6  nm  were  used.30 
The  error  in  the  temperature  determined  by  utilizing  these  two  lines  is  small,  less  than  5%  for  the 
entire  range  from  500K  to  2500K.31  Temperatures  were  derived  by  comparing  the  relative  integrated 
intensity  of  these  two  lines  in  the  experimental  spectra  to  those  in  spectra  predicted  at  various  temper¬ 
atures  using  the  computer  program  LIFBASE.32  The  pressure  was  put  into  the  LIFBASE  program 
and  Doppler,  predissociation,  and  collisional  broadening  were  included  in  the  spectral  simulations. 
The  instrument  resolution  parameter  was  adjusted  until  the  predicted  peak  width  matched  our 
observed  data,  which  occurred  at  an  instrument  resolution  of  0.035 A.  Then  a  series  of  simulated 
spectra  were  run  for  temperatures  over  the  expected  experimental  temperature  range  of  500-2000K. 

A  relation  for  temperature  as  a  function  of  the  ratio  of  the  areas  under  the  Pi(2)  and  Rj(l  4)  peaks  was 
obtained.  This  function  was  then  used  to  generate  a  temperature  for  each  experimental  spectrum  by 
using  the  integrated  intensities  of  the  observed  Pi(2)  and  Ri(14)  peaks  after  correcting  for  laser 
power.  The  relative  OH  concentration  was  measured  by  comparing  the  integrated  intensity  of  each  of 
these  lines  after  correcting  for  the  temperature  dependence.  This  method  gave  OH  concentration  pro¬ 
files  similar  to  those  obtained  by  monitoring  the  Pi(8)  line  at  285.66  nm,  which  has  been  shown  to  be 
relatively  insensitive  to  temperature  from  1200K  to  2500K.9 

In  these  experiments,  we  measured  relative  OH  concentrations.  Absolute  concentrations  were  not 
measured  due  to  the  complexity  of  accounting  for  excited-state  quenching.  In  order  to  obtain  an 
absolute  concentration  without  making  a  quantitative  absorption  measurement  for  calibration,  the 
identity,  concentration,  and  quenching  cross  section  of  each  species  that  might  collide  with  the  OH 
must  be  known  at  every  point.  In  the  flames  with  no  diluent,  the  main  quenching  species  is  water,  but 
in  the  flames  with  argon  diluent,  the  argon  also  plays  an  important  role  in  quenching.  Argon  has  a 
quenching  cross  section  for  OH  that  is  several  orders  of  magnitude  smaller  than  that  of  water.33 
Because  the  quenching  environment  is  so  different  in  the  experiments  with  and  without  diluent,  the 
relative  OH  signals  for  experiments  with  diluent  cannot  be  directly  compared  to  those  in  experiments 
without  diluent.  The  addition  of  argon  increased  the  width  of  the  time  evolution  of  the  OH  fluores¬ 
cence  signal,  especially  in  the  tail.  Based  on  the  changes  in  the  width  and  shape  of  this  signal  as 
argon  was  added,  it  seems  to  be  very  sensitive  to  the  quenching  rate.  The  addition  of  small  amounts 
of  SiCl4  did  not  perturb  the  width  or  shape  of  the  time  evolution  of  the  OH  fluorescence  signal.  Thus, 
we  assume  that  the  change  in  the  quenching  environment  on  addition  of  small  amounts  of  SiCl4  is 
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sufficiently  small  that  we  can  compare  OH  concentrations  for  experiments  with  the  same  basic  flame 
configuration  with  and  without  SiCU  precursor. 

Hydroxyl  concentration  and  temperature  were  monitored  as  a  function  of  distance  from  the  burner.  A 
burner  position  of  “0.0  cm”  indicates  that  the  burner  blocks  half  of  the  0.6-mm-dia  laser  beam.  The 
burner  was  moved  through  several  centimeters  for  each  data  run,  with  data  being  collected  at  various 
distances  from  the  burner.  Although  we  collected  data  out  to  6  cm  from  the  burner  for  some  runs, 
most  runs  ended  at  3  cm,  at  which  point  combustion  is  substantially  complete.  Beyond  3  cm  from  the 
burner,  the  temperature  and  OH  concentration  depend  mostly  on  the  rate  of  heat  loss  from  the 
apparatus,  which  does  not  shed  any  light  on  the  chemistry.  The  error  bars  shown  in  the  figures 
represent  one  standard  deviation  based  on  several  data  runs. 
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3.  Mechanism  and  Model 


Table  1  outlines  our  proposed  combustion  mechanism  for  SiCI4  in  a  hydrogen/oxygen  flame.  Rates 
for  the  first  thermal  decomposition  step  and  the  first  hydrogen  abstraction  step  were  available  in 

Table  1.  Silicon  Tetrachloride  Proposed  Combustion  Mechanism  (k  =  ATbe  BRT).  The  units 
of  A  are  cc/molecule-s  for  second-order  reactions;  higher  order  reactions  have 
higher  power  inverse  concentration  units.  The  activation  energy  E  is  in  kJ/mol. 


Reaction 

A 

b 

E 

Source 

SiCL,  +  M  =  SiCI3  +  Cl  +  M 

4.8e16 

0 

392 

Experiment34 

SiCI,  +  H  =  SICI3  +  HCI 

1.7e13 

0 

29 

This  work 

Original  rate: 

1.4e13 

0 

40 

Experiment35 

SiCI3  +  O2  +  M  =  SiCI302  +  M 

8.7e12 

0 

0 

Experiment1361 

SiCI302  +  M  ss  SiCI30  +  O  +  M 

1.0e13 

0 

59 

Analogy151 

SiCI30  +  H  =  SiCI20  +  HCI 

1.6e17 

-1.4 

6 

Analogy42 

SiCI302  +  H02  ss  SiCI302H  +  02 

4.0e9 

0 

0 

Analogy151 

SiCI20  +  M  =  SiO  +  Cl2  +  M 

1.0e12 

0 

254 

Analogy42 

SiCIjO  +  H  =  SiCIO  +  HCI 

2.0e14 

0 

21 

This  work 

Original  rate: 

3.3e13 

0 

44 

Analogy151 

SICIO  +  M  =  SiO  +  Cl  +  M 

1.0el4 

0 

33 

This  work 

Original  rate: 

5.4e13 

0 

121 

Analogy151 

SiCIO  +  H  =  SiO  +  HCI 

3.0e15 

0 

0 

This  work 

Original  rate: 

2.0e13 

0 

0 

Analogy151 

SiCI3  +  H  =  SiCI2  +  HCI 

1.5e12 

0 

10 

Analogy141 

SiCls  +  H  =  SiCI  +  HCI 

1.4e12 

0 

8 

Analogy42 

SiCI2  +  02  =  SiCI20  +  0 

3.4e12 

0 

0 

Experiment1371 

SiCI2  +  0  ~  SiO  +  Cl2 

4.6e12 

0 

0 

Analogy1381 

SiCI +  0  =  SiO  +  Cl 

4.0e12 

0 

0 

Analogy42 

SiCI30  +  H2  =  SiCI3OH  +  H 

1.2e11 

0 

4 

Analogy151 

SiCI302  +  H2  =  SiCI302H  +  H 

1.3e10 

0 

36 

Analogy151 

SiCI302H  +  M  =  SiCI30  +  OH  +  M 

6.5e13 

0 

204 

Analogy151 

SiCI3OH  +  OH  =  SiCI30  +  H20 

4.0e11 

0 

6 

Analogy151 

SiCI20  +  H20  =  SiC!02H  +  HCI 

4.0e12 

0 

81 

Calculated1’11 

SiCI02H  +  M  =  Si02  +  HCI  +  M 

3.7e13 

0 

205 

Calculated1’’1 

SiCI20  +  OH  =  SiCI02H  +  Cl 

1.0e12 

0 

0 

Analogy151 

SiCI3  +  SiOI3  =  SiCl2  +  SiCI4 

1.0e12 

0 

0 

Fit1391 

SiCI3  +  M  =  SiCI2  +  Cl  +  M 

1.0e14 

0 

272 

Fit1391 

SiCI4  +  H20  =  SiCI3OH  +  HCI 

1.0e12 

0 

121 

Experiment1401 

SiCI3OH  +  M  =  SiCI20  +  HCI  +  M 

2.8e13 

0 

207 

Calculated1’’1 

SiCI4  +  OH  =  SiCI3OH  +  Cl 

1.4e7 

0 

188 

Calculated1”1 

S10  +  02  =  Si02+0 

l.5e12 

0 

21 

This  work 

Reaction 

A 

b 

E 

Source 

Original  rate: 

1.0e13 

0 

27 

Estimated14’1 

SiO  +  OH  =  SI02  +  H 

3.0et0 

0.8 

0 

This  work 

Original  rate: 

1.8e10 

0.78 

5 

Calculated'31 

SiO  +  0  +  M  =  Si02  +  M 

2.5e15 

0 

18 

Analogy  with  CO1411 

HCI  +  M  =  H  +  Cl  +  M 

4.4e13 

0 

342.1 

NIST  database43 

Cl2  +  M  =  2CI  +  M 

2.32e13 

0 

196.5 

NIST  database43 

H  +  Cl2  =  HCI  +  Cl 

1.91e13 

0 

0 

NIST  database43 

0  +  HCI  =  OH  +  Cl 

6.03e12 

1.7 

27.4 

NIST  database43 

OH  +  HCI  =  H20  +  Cl 

1.11e7 

0 

-2.8 

NIST  database43 

papers  by  Catoire  et  a/.34'35  and  Powers,14  respectively.  These  papers  suggested  that  the  initial  steps 
produce  the  radical  SiCl3  rather  than  SiCl2.  Because  the  analogous  radical  SiH3  is  an  important  spe¬ 
cies  in  the  Takahashi  et  al.  mechanism5  for  silane  oxidation,  we  used  this  mechanism  as  a  model  for 
our  silicon  tetrachloride  oxidation  mechanism.  Additional  reactions  for  species  in  the  mechanism 
were  added  whenever  experimental,36’37'38’39’40  calculated,  or  estimated41  rates  were  available.  The 
reactions  were  chosen  as  analogues  of  the  corresponding  reactions  for  silane,  with  the  exception  that 
no  chlorine  oxides  were  included  since  these  are  not  stable  enough  to  exist  at  the  flame  temperatures 
of  interest.  Also,  we  did  not  include  any  species  in  which  hydrogen  re-substituted  on  the  silicon  since 
this  is  unfavorable  energetically.  Analogous  rates  were  taken  from  the  Takahashi  list5  or  the  Zacha- 
riah  mechanism.42  When  estimating  chlorine  rates  from  analogous  hydrogen  rates,  rates  were  esti¬ 
mated  as  roughly  10  times  slower  for  abstraction,  dissociation,  and  transfer  reactions,  and  roughly  the 
same  for  decomposition  and  reactions,  which  required  more  indirect  analogies.  These  estimates  are 
based  on  the  few  reactions  where  experimental  rate  determinations  have  been  made  for  both  the  silane 
rate  and  the  analogous  silicon  chloride  rate.  There  is  not  enough  information  to  justify  a  more  precise 
estimate.  Because  this  estimate  is  so  crude,  we  tried  adjusting  it  by  a  factor  of  10  in  either  direction 
during  the  fitting  process  to  try  to  get  better  agreement  with  the  experiment.  Adjusting  the  rates  as  a 
group  did  not  change  the  fit,  but  changing  some  of  the  individual  rates  in  this  group  was  helpful.  HC1 
and  CI2  reactions  were  added  as  necessary  from  the  NIST  database.43  The  GRI  3.0  mechanism44  was 
used  for  the  H2/02  reactions.  Particle  formation  and  chemistry  on  particle  surfaces  were  not  included 
in  the  model.  Because  particle  formation  is  a  3 -body  process,  it  is  expected  to  occur  relatively  late  in 
the  flame.  Because  of  diffusion,  however,  particles  are  distributed  throughout  the  flame,  and  there  is 
some  potential  for  reactions  on  the  particle  surfaces  to  influence  flame  chemistry. 

In  three-body  reactions,  the  reaction  is  dependant  on  a  non-reacting  collider  molecule,  or  “third 
body,”  to  transfer  energy  to  or  from  the  reacting  species.  Because  energy  transfer  and  storage  vary 
greatly  for  different  chemical  species,  different  third-body  colliders  have  different  efficiencies  at 
promoting  the  reaction.  The  GRI  mechanism  accounts  for  the  difference  in  efficiency  of  various  third 
bodies;  for  instance,  argon  is  usually  1-3  orders  of  magnitude  less  efficient  than  water  as  a  third  body. 
Based  on  this,  we  have  also  reduced  the  efficiency  of  argon  as  a  third  body  in  the  silicon  reactions  by 
2  orders  of  magnitude. 

We  compared  our  data  to  results  from  a  chemical  kinetics  model  using  our  proposed  SiCU  combus¬ 
tion  mechanism  and  CHEMKIN.  We  chose  the  premixed  burner  reactor  in  the  CHEMKIN  program 
to  model  our  one-dimensional  flame.  This  model  uses  the  measured  temperature  profile  as  an  input. 
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We  did  a  fit  to  our  temperature  data  and  used  this  fit  as  the  model  input.  For  each  set  of  basic  flame 
configurations,  we  scaled  the  measured  relative  OH  concentrations  to  the  model  prediction  for  the  OH 
concentration  of  the  H2/O2  flame  with  no  precursor. 

The  CHEMKIN  transport  subroutine  requires  input  parameters  to  calculate  the  diffusivity  of  the 
various  species.  We  estimated  these  parameters  for  silicon  tetrachloride  combustion  species  based  on 
the  methods  outlined  in  “The  Properties  of  Gases  and  Liquids”  by  Reid,  Prausnitz,  and  Sherwood.45 
Values  for  the  dipole  moments  came  from  the  CRC  handbook46  where  available;  calculated  values 
from  Gaussian47  were  used  where  no  measured  values  were  available.  The  transport  properties  are 
shown  in  Table  2. 

Table  2.  Transport  Database  for  Species  from  SiCl4  Combustion.  ()  =  estimated  quantity 


Species  Name 

Geometry# 

e/k 

a 

P 

HSiO 

1 

500 

3.69 

(2) 

SiOH 

1 

363 

4.11 

(1) 

SiH20 

2 

500 

3.71 

(2) 

SiO 

1 

3150 

2.96 

3.1 

Si02 

1 

4300 

3.51 

0 

H3SiOO 

2 

380 

4.54 

(1.5) 

HsSiOO 

2 

380 

4.51 

(1.5) 

OSiOH 

1 

435 

4.28 

(2) 

HjSiO 

2 

4400 

1.80 

4.6 

H3SiO 

2 

2600 

2.17 

(4) 

H3SiOH 

2 

400 

4.05 

(1.5) 

HjSiOH 

2 

480 

3.77 

(2) 

H3SiOOH 

2 

380 

4.56 

(1.5) 

OSiHOH 

2 

350 

4.61 

(0) 

H3SiOO* 

2 

380 

4.54 

(1.5) 

SiOOH 

2 

350 

4.59 

(0) 

SiH„ 

2 

190 

4.21 

0 

HSiOH 

1 

363 

4.14 

(1) 

HSiOOH 

1 

435 

4.31 

(2) 

Si202 

2 

2360 

4.43 

(0) 

Si2H 

2 

236 

4.49 

(0.1) 

H3Si2H 

2 

300 

5.19 

(0.1) 

SiCI 

1 

366 

4.01 

1.03 

SiCI2 

1 

360 

4.69 

0.95 

SiCI3 

2 

350 

5.22 

0.74 

SiCI4 

2 

390 

5.64 

0 

SiCI302 

2 

473 

5.60 

(1) 

SiCI30 

2 

473 

5.41 

(1) 

SiCI20 

2 

970 

3.90 

3.79 

SiCI302H 

2 

473 

5.62 

(1) 

SiCI3OH 

2 

473 

5.42 

(1) 

SiCI02H 

2 

476 

4.65 

(1) 

SiCIO 

1 

556 

4.15 

(2) 

#:  Geometry  index.  0  =  monatomic,  1  =  linear,  2  =  nonlinear 
e/k,  Lennard-Jones  well  depth  in  Kelvins 
a,  Lennard-Jones  collision  diameter  in  Angstroms 
p,  Dipole  moment  in  Debye  units. 
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The  modeling  work  required  additional  thermodynamic  data  for  the  CHEMKIN  database.  The  ther¬ 
modynamic  data  are  necessary  to  calculate  the  reverse  rate  of  reactions.  This  database  consists  of  sets 
of  coefficients  that  represent  the  entropy,  enthalpy,  and  heat  capacity  of  each  substance  as  a  power 
series  expansion  in  terms  of  temperature.  The  expansion  used  in  CHEMKIN  is  the  Gordon  and 
McBride  representation.24  Most  of  the  literature  values,  including  the  NIST  database,43  are  given  as 
Shomate  coefficients,  which  are  for  a  similar,  but  not  identical,  expansion.  We  were  able  to  relate  the 
two  sets  of  coefficients  by  approximation.  Enthalpies  as  a  function  of  temperature  calculated  from 
the  converted  coefficients  were  compared  with  the  JANAF  Tables48  where  possible.  This  approach 
was  found  to  give  results  consistent  to  about  1%  for  molecules  in  the  NIST  database.  For  species  not 
in  the  NIST  database,  the  estimated  thermodynamic  quantities  were  taken  from  formal  calculations, 
such  as  those  in  Zachariah  and  Tsang3  or  Allendorf,  et  al.49  For  species  not  in  the  NIST  database  and 
for  which  no  formal  calculations  exist,  we  use  the  three-term  representation: 

ai  =  Cp/R 

a6  =  H°/R  -  a)(298) 

a7  =  S0/R-a,ln(298) 

In  addition,  we  used  approximations  for  the  heat  capacity  and  the  entropy  terms  for  molecules  where 
these  are  not  available.  For  the  heat  capacity,  we  use  the  classical  mechanics  approximation:  Cp  = 
(3n  -  2)R,  where  n  is  the  number  of  atoms  in  the  molecule.  For  the  entropy,  we  use  an  empirical 
correlation  with  the  molecular  weight:  S°  =  167  +  2.163(M.W.)  J/mol/K.  The  enthalpy  was  based 
on  approximate  bond-additivity  relations.  The  thermodynamic  coefficients  are  listed  in  Table  3. 


Table  3.  Thermodynamic  Data  for  the  Silicon  Species  in  Mechanism.  Species  with  a5  =  0  are 
from  NIST  database.  Species  with  all  coefficients  taken  from  calculations.  Other 
species  coefficients  estimated  from  heat  capacity,  enthalpy,  and  entropy  at  298K. 


Species 

a, 

fl2 

a3 

a. 

a5 

a8 

a7 

Si02(g) 

7.4037 

5.37E-05 

-1.09E-08 

7.52E-13 

0.0 

-4.01  E+04 

-1.63E+01 

Si02(l) 

10.316 

-1.92E-9 

4.81E-13 

-4.58E-17 

0.0 

-1.15E+05 

-5.76E+01 

SiO(g) 
low  T 

2.3482 

4.51  E-03 

-3.67E-06 

1.09E-09 

0.0 

-1.29E+04 

1.10E+01 

SiO(g) 

highT 

4.2936 

2.08E-04 

-6.13E-08 

7.14E-12 

0.0 

-1.38E+04 

3.12E-01 

81(9) 

1.7553 

6.28E-04 

-1 .30E-07 

8.90E-12 

0.0 

5.42E+04 

1.08E+01 

SiCI302 

1.61  E+01 

0 

0 

0 

0 

-7.58E+04 

-4.14E+01 

SiCI302H 

1.91E+01 

0 

0 

0 

0 

-9.68E+04 

-5.61  E+01 

SiCI3OH 

highT 

1.27E+01 

3.13E-03 

-1.33E-06 

2.80E-10 

-2.35E-14 

-1.05E+05 

-3.10E+01 

SiCI3OH 

lowT 

4.77E+00 

3.86E-02 

-6.25E-05 

4.79E-08 

-1.40E-1 1 

-1.03E+05 

6.31  E+00 

SiCI30 

highT 

1.06E+01 

2.06E-03 

-1.21E-06 

3.29E-10 

-3.40E-14 

-6.64E+04 

-2.04E+01 

SiCI30 

lowT 

4.27E+00 

2.98E-02 

-4.87E+03 

3.74E-08 

-1.10E-1 1 

-6.52E+04 

9.53E+00 

SiCI20 
high  T 

8.25E+00 

2.52E-03 

-1.49E-06 

4.10E-10 

-4.28E-14 

-5.87E+04 

-1.17E+01 
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Species 

a, 

a2 

a3 

as 

as 

a7 

SiCI20 

lowT 

4.27E+00 

1.90E-02 

-2.84E-05 

2.07E-08 

-5.92E-12 

-5.79E+04 

7.42E+00 

SiCIO 

highT 

5.72E+00 

1.86E-03 

-1.12E-06 

3.10E-10 

-3.27E-14 

-2.13E+04 

1.00E-01 

SiCIO 

lowT 

3.77E+00 

9.20E-03 

-1.20E-05 

7.80E-09 

-2.06E-12 

-2.08E+04 

9.66E+00 

SiCI02H 

highT 

8.76E+00 

4.49E-03 

-2.16E-06 

5.15E-10 

-4.87E-14 

-8.01  E+04 

-1.56E+01 

SiCI02H 

lowT 

2.52E+00 

3.08E-02 

-4.51  E-05 

3.24E-08 

-9.10E-12 

-7.89E+04 

1.43E+01 

SiCI302H 

1.91E+01 

0 

0 

0 

0 

-9.68E+04 

-5.61  E+01 
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4.  Results  and  Analysis 


Figure  2  shows  typical  fluorescence  data  for  five  data  runs  at  a  distance  of  1.2  cm  from  the  burner. 
These  data  have  been  corrected  for  laser  power.  The  symbols  represent  the  experimental  data,  and 
the  lines  are  the  LIFB  ASE  simulated  spectra  (normalized  to  the  data)  for  the  temperatures  given  by 
the  relative  intensities  of  the  experimental  Pi(2)  and  Ri(14)  peaks.  Below  about  900K,  the  signal-to- 
noise  ratio  for  the  Ri(14)  peak  is  less  than  one,  and  although  the  integrated  intensity  can  still  be 
obtained,  the  error  bars  are  higher. 


Figures  3  and  4  show  temperature  and  OH  concentration,  respectively,  for  three  basic  flame  configu¬ 
rations  with  and  without  SiCl4:  an  undiluted  stoichiometric  H2/O2  flame,  a  stoichiometric  H2/O2 
flame  diluted  with  argon,  and  an  undiluted  rich  H2/02  flame  (<b  =  1.2).  The  error  bars  in  these  figures 
are  likely  due  mainly  to  a  combination  of  shot-to-shot  instabilities  inherent  in  a  Nd:YAG  laser  and  the 
slight  difference  in  the  environment  for  each  shot  due  to  particles. 


In  general,  the  flame  temperature  increased  upon  addition  of  SiCl4.  This  effect  was  the  largest  for  the 
stoichiometric  undiluted  flame,  as  seen  in  Figure  3A,  where  the  maximum  flame  temperature  without 
SiCU  was  1 125K  but  increased  by  180K  to  1300K  when  SiCU  was  added.  For  the  diluted  flame 
(Figure  3B),  the  temperatures  were  slightly  lower  and  rose  more  slowly  up  to  a  maximum  of  1 120K 
without  SiCU  and  1200K  with  SiCU,  a  temperature  rise  of  80K.  For  the  rich  flame  (Figure  3C),  the 
temperature  difference  upon  adding  SiCU  is  even  less,  reaching  1060K  for  the  flame  without  SiCl4 
and  only  increasing  35K  to  1095K  for  the  flame  with  SiCl4. 


wavelength,  nm 

Figure  2.  Fluorescence  data  collected  1.2  cm  from  the  burner  for  five  data  runs.  The 

symbols  represent  the  experimental  data,  and  the  lines  represent  the  LIFBASE 
simulated  spectrum  for  the  temperature  given  by  the  ratio  of  the  areas  under  the 
P,(2)  and  Ri(14)  peaks.  The  circles  and  solid  line  are  for  run  1  (temperature 
129 IK),  the  triangles  pointing  down  and  dotted  line  are  for  run  2  (temperature 
1307K),  the  squares  and  short-dashed  line  are  for  run  3  (temperature  1253K),  the 
diamonds  and  dot-dot-dashed  line  are  for  run  4  (temperature  1 320K),  and  the 
triangles  pointing  up  and  the  long-dashed  line  are  for  run  5  (temperature  1260K). 
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Distance  from  burner  face,  mm 

Figure  3.  Temperature  profiles  for  hydrogen/oxygen  flames  with  and  without  SiCl4.  (A) 

shows  data  for  the  stoichiometric  undiluted  flame.  The  circles  indicate  data  with¬ 
out  SiCl4  and  the  squares  show  data  for  the  flame  with  SiCl4.  (B)  shows  data  for 
the  stoichiometric  diluted  flame.  The  diamonds  indicate  data  without  SiCl4  and 
the  hexagons  show  data  for  the  flame  with  SiCl4.  C)  shows  data  for  the  rich  undi¬ 
luted  flame.  The  triangles  pointing  up  indicate  data  without  SiCl4  and  the  triangles 
pointing  down  show  data  for  the  flame  with  SiCl4.  The  lines  indicate  fits  to  the 
temperature  data. 
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Distance,  mm 

Figure  4.  OH  concentration  profiles  for  hydrogen/oxygen  flames  with  and  without  SiCl4.  In 
all  panels,  the  solid  lines  are  the  model  predictions  of  OH  concentration  for  the 
flames  without  SiCl4,  the  dashed  lines  are  the  model  predictions  of  the  OH  con¬ 
centration  for  the  flames  with  SiCl4  using  the  modified  rates,  and  the  dotted  lines 
are  the  model  predictions  of  the  OH  concentration  for  the  flames  with  SiCl4  using 
the  original  rates.  A)  shows  data  for  the  stoichiometric  undiluted  flame.  The  cir¬ 
cles  indicate  data  without  SiCl4  and  the  squares  show  data  for  the  flame  with  SiCl4. 

B)  shows  data  for  the  stoichiometric  diluted  flame.  The  diamonds  indicate  data 
without  SiCl4  and  the  hexagons  show  data  for  the  flame  with  SiCl4.  C)  shows  data 
for  the  rich  undiluted  flame.  The  triangles  pointing  up  indicate  data  without  SiCl4 
and  the  triangles  pointing  down  show  data  for  the  flame  with  SiCl4. 

The  addition  of  SiCU  suppressed  the  OH  concentration  for  all  flame  configurations,  as  can  be  seen  in 
Figure  4.  For  the  undiluted  flame  (Figure  4A),  the  addition  of  SiCl4  reduced  the  OH  concentration  to 
-85%  of  its  value  without  SiCU.  In  the  diluted  flame  (Figure  4B),  the  reduction  in  OH  concentration 
on  addition  of  SiCU  was  similar,  with  the  SiCU  flame  having  -79%  of  the  OH  of  the  flame  without 
SiCU-  However,  in  the  rich  flame  (Figure  4C)  the  reduction  was  much  higher,  with  the  SiCU  flame 
having  only  35%  of  the  maximum  OH  observed  in  the  flame  without  SiCU-  The  shape  of  the  OH 
concentration  profile  changed  more  on  adding  SiCU  for  the  rich  flame,  becoming  much  more  flat  than 
the  OH  profile  without  SiCl4. 
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The  model  prediction  of  the  OH  concentration  for  the  stoichiometric  undiluted  flame  without  SiCl4 
was  within  the  error  bars  of  the  measured  values  except  at  5-7  mm  above  the  burner,  where  the  model 
slightly  underpredicted  the  OH  concentration.  The  fit  to  the  rich  flame  without  SiCl4  was  good  for 
the  leading  edge  and  tail  of  the  OH  concentration  peak,  but  OH  was  underpredicted  in  the  region  from 
15  to  20  mm  above  the  burner.  The  model  prediction  of  the  OH  concentration  for  the  diluted  flame 
without  SiCl4  significantly  underpredicted  the  OH  concentration  for  the  entire  leading  edge  of  the  OH 
peak,  the  region  from  the  burner  surface  to  12  mm  above  the  burner.  Even  adjusting  the  input  tem¬ 
perature  profile  to  the  extremes  of  the  error  bars  did  not  correct  this  problem.  All  three  of  these  H2/O2 
flames  without  SiCl4  were  modeled  using  the  GRI  3.0  mechanism  and  the  premixed  burner  input  to 
CHEMKIN,  so  the  comparison  of  the  prediction  to  the  experimental  data  gives  some  idea  of  the 
limitations  of  this  model. 

The  model  prediction  of  OH  concentration  for  the  stoichiometric  hot  flame  with  SiCl4  fits  well  in  the 
leading  edge  and  peak  of  the  OH  concentration,  but  underpredicts  the  OH  in  the  tail,  the  region  more 
than  14  mm  from  the  burner.  For  the  rich  flame  without  precursor,  the  shape  of  the  predicted  OH 
profile  is  somewhat  off,  overpredicting  OH  in  the  early  part  of  the  flame  from  4  to  1 0  mm  from  the 
burner  and  underpredicting  in  the  tail  (beyond  15  mm).  The  experimental  OH  profile  is  more  flat 
than  the  model  prediction.  However,  the  magnitude  of  the  OH  concentration  is  correctly  predicted. 
The  model  significantly  underpredicts  the  OH  concentration  throughout  the  entire  region  for  the 
diluted  flame  with  precursor. 
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5.  Discussion 


5.1  Reaction  Rates 

The  initial  set  of  rates  chosen  for  our  investigation  of  the  SiCl4  mechanism  did  not  provide  a  satis¬ 
factory  fit  to  the  data  for  any  of  the  flame  configurations  studied,  as  can  be  seen  by  the  dotted  lines  in 
Figure  4.  It  overpredicted  the  OH  concentration  by  a  factor  of  3  for  the  undiluted  stoichiometric 
flame  (Figure  4A),  a  factor  of  2  for  the  diluted  flame  (Figure  4B),  and  a  factor  of  4.5  for  the  rich 
flame  (Figure  4C).  We  performed  a  sensitivity  analysis  in  CHEMKIN  to  determine  which  rates 
affected  the  OH  profile  the  most.  The  sensitivities  of  the  OH  concentration  to  various  rates  for  the 
final  rate  set  are  shown  in  Table  4.  While  several  different  combinations  of  rates  could  yield  a  rea¬ 
sonable  prediction  of  the  OH  concentration  for  the  stoichiometric  undiluted  flame,  we  chose  the  com¬ 
bination  that  also  did  the  best  job  at  fitting  the  data  from  the  other  two  flame  configurations.  The 
predicted  OH  concentration  was  sensitive  to  several  reactions  that  do  not  have  well-known  rates. 
Table  5  lists  these  reactions  and  the  modifications  we  made  to  the  rate  expressions. 

Table  4.  CHEMKIN  Sensitivities  of  OH  Concentration  to  the  Reaction  Rates  at 
Different  Distances  from  the  Burner.  The  sensitivities  are  listed  as  fol¬ 


lows:  the  undiluted  stoichiometric  flame  in  bold,  the  diluted  stoichio¬ 
metric  flame  in  plain  text,  and  the  undiluted  rich  flame  in  italics. 


Reaction 

Distance  from  burner  (cm) 

0.5  1.0  2.0  3.0 

OH  +  H2  =  H  +  H20 

-0.797 

-0.520 

-0.328 

-0.288 

-0.940 

-0.844 

-0.706 

-0.695 

-0.948 

-0.944 

-0.959 

-0.957 

H  +  02  =  0  +  OH 

0.349 

0.115 

-0.079 

-0.080 

0.135 

0.262 

0.046 

-0.037 

0.151 

0.034 

-0.134 

-0.175 

Cl  +  H2  =  HCI  +  H 

0.060 

0.014 

0.008 

0.012 

-0.087 

0.043 

0.000 

-0.004 

0.045 

0.006 

-0.015 

-0.022 

20H  =  0  +  H20 

-0.018 

-0.035 

-0.051 

-0.056 

-0.009 

-0.013 

-0.018 

-0.015 

-0.002 

-0.001 

0.000 

0.001 

H  +  H02  =  20H 

0.010 

0.000 

-0.002 

-0.002 

0.031 

0.006 

0.000 

-0.002 

0.012 

-0.001 

-0.011 

-0.014 

SiO  +  02  =  Si02  +  0 

0.484 

0.260 

0.046 

-0.102 

0.698 

0.596 

0.469 

0.375 

0.672 

0.316 

-0.478 

-1.325 

SiCIO  +  M  =  SiO  +  Cl  +  M 

-0.230 

-0.130 

0.016 

0.033 

-0.291 

-0.198 

-0.072 

-0.010 

-0.168 

-0.041 

0.132 

0.214 

SiCIO  +  H  =  SiO  +  HCI 

-0.095 

-0.114 

-0.109 

-0.104 

-0.172 

-0.092 

-0.108 

-0.114 

-0.121 

-0.045 

0.069 

0.173 
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Reaction 

Distance  from  burner  (cm) 

0.5  1.0  2.0  3.0 

SiCI4  +  H  =  SiCI3  +  HCI 

-0.014 

-0.024 

0.014 

0.019 

0.012 

0.017 

0.091 

0.039 

-0.063 

-0.184 

SiO  +  OH  =  Si02  +  H 

0.005 

■ns 

0.000 

0.347 

0.012 

0.012 

0.078 

-0.017 

-0.078 

-0.139 

SiCI20  +  H  =  SiCIO  +  HCI 

0.007 

1TT1 

0.036 

0.018 

0.024 

0.028 

0.029 

0.015 

-0.016 

-0.068 

SiCI3  +  H  =  SiCI2  +  HCI 

-0.003 

-0.013 

0.012 

0.022 

0.011 

-0.018 

-0.048 

SiCI20  +  H20  =  SiCI02H  +  HCI 

-0.004 

mmm 

-0.010 

-0.016 

-0.019 

-0.016 

-0.013 

0.048 

SiCI20  +  OH  =  SICI02H  +  Cl 

-0.003 

-0.012 

0.015 

Table  5.  Reactions  modified  in  order  to  improve  the  model  fit  to  our  data 


Reaction 

Source 

Modification 

SiO  +  02  =  Si02  +  O 

Estimated  by  Jachimowski41 

Decrease  A  by  ~  order  of  magnitude 
Decrease  Ea  from  27  to  21  kJ/mol 

SiCIO  +  H  =  SiO  +  HCI 

Estimated  by  analogy  with  SiH4 15! 

Increase  A  by  factor  of  150 

Ea  was  estimated  to  be  0;  unchanged 

SiO  +  OH  =  SiO,  +  H 

Calculated  by  Zachariah131 

Increase  A  by  ~  factor  of  2 

Decrease  Ea  from  5  to  0  kJ/mol 

SiCIO  +  M  =  SiO  +  Cl 

Estimated  by  analogy  with  SiH4 151 

Increase  A  by  factor  of  2 

Decrease  Ea  from  121  to  33  kJ/mol 

SiCI4  +  H  =  SiCI3  +  HCI 

Shock  tube  experiment  by  Catoire  et  al.  34 

Increase  A  by  ~  25% 

Decrease  Ea  from  40  to  29  kJ/mol 

SiCI20  +  H  =  SiCIO  +  HCI 

Estimated  by  analogy  with  SiH„ ,5) 

Increase  A  by  ~  order  of  magnitude 
Decrease  Ea  from  44  to  21  kJ/mol 

The  rate  for  the  reaction  to  which  the  predicted  OH  concentration  profile  was  most  sensitive,  SiO  + 
02,  was  thought  to  be  too  high  by  Babushok  et  al.,  who  estimated  an  activation  energy  (Ea)  of  125 
kJ/mol  and  an  Arrhenius  parameter  (A)  of  5  x  1012  cc/molecule-s.  When  we  used  these  parameters 
and  adjusted  the  other  reaction  parameters  until  the  predicted  peak  had  the  observed  magnitude  for 
the  stoichiometric  undiluted  flame,  the  OH  concentration  was  significantly  underpredicted  in  the  ris¬ 
ing  edge  and  in  the  tail.  The  predictions  for  the  other  two  flame  configurations  were  also  poor.  The 
original  estimate  of  Ea  =  27  kJ/mol  and  A  =1  X  1013  cc/molecule-s  was  made  by  Jachimowski,  who 
used  Hirschfelder’s  rule  to  estimate  the  activation  energy  and  analogy  to  CO  to  determine  the  A  fac¬ 
tor.  We  reduced  the  activation  energy  and  the  Arrhenius  parameter  from  Jachimowski’s  estimate  to 
fit  our  data. 

Of  the  six  reactions  that  we  modified  to  fit  our  data,  one  has  been  experimentally  measured.  Catoire 
et  al?5  measured  the  reaction  rate  of  SiCl4  +  H  in  a  shock  tube  experiment  at  1530  and  1730K  and 
obtained  an  Arrhenius  parameter  of  A  =  1.4  x  1013  (±24%)  cc/molecule-s  and  an  activation  energy  of 


18 


Ea  =  40±22  kJ/mol.  Our  values  for  the  activation  energy  and  Arrhenius  parameter  fall  within  the 
reported  experimental  error  bar  range.  This  rate  was  also  calculated  by  Zhang  et  al.,50  who  calculated 
A  =  6.6  x  1014  cc/molecule-s  and  Ea=  1 13  kJ/mol  for  the  temperature  range  1 500-1 800K.  The  rate 
calculated  by  Zhang  is  within  the  error  of  the  experimental  measurement  over  the  temperature  range 
in  which  the  measurements  were  made.  Some  reasons  for  the  large  difference  in  activation  energy 
between  the  experiment  and  the  calculation  are  discussed  in  Zhang’s  paper50  Another  consideration  is 
that  the  basis  set  used  in  the  calculations,  QCISD,  has  been  found  to  give  higher  reaction  energies  and 
barrier  heights  than  other  basis  sets  for  reaction  of  chlorine  with  methane.51  The  value  for  this  rate 
that  we  used  in  our  model  in  order  to  provide  the  best  fit  to  our  data  is  closer  to  the  experimental  rate 
than  to  the  calculated  one. 


5.2  Temperature 

Addition  of  SiCl4  precursor  caused  the  flame  temperature  to  increase  for  all  three  basic  flame  con¬ 
figurations.  An  equilibrium  analysis  based  on  our  model  indicates  that  there  is  very  little  energy  in 
the  combustion  of  SiCl4  compared  to  H2  and  O2.  In  fact,  combining  with  H2  and  O2  to  form  HC1  and 
SiCbCg)  releases  significantly  less  heat  than  if  the  H2  and  O2  reacted  to  form  H2O.  The  only  signifi¬ 
cant  heat  release  comes  from  the  condensation  of  the  silicon  oxides;  this  brings  the  total  heat  release 
up  to  about  the  same  amount  as  for  water  formation.  The  heat  release  from  particle  formation  is 
noticeable  because  it  occurs  in  a  portion  of  the  flame  different  from  the  heat  release  from  water  for¬ 
mation.  Therefore,  the  region  of  the  flame  where  the  temperature  increase  is  observed  could  indicate 
the  region  of  particle  formation.  SiCl4  has  a  higher  heat  capacity  than  the  diatomics  H2  and  O2;  this 
causes  a  slower  rise  in  temperature  in  the  early  portion  of  the  SiCl4-doped  flame. 

A  flame  temperature  rise  upon  precursor  addition  was  also  observed  by  Hwang  et  al .7  as  well  as  by 
Wilcox  et  al ,52  Hwang  et  al.  attributed  the  increase  of  ~100K  in  the  region  35  to  80  mm  above  the 
burner  to  heat  released  during  exothermic  reactions  of  SiCl4  oxidation  and  hydrolysis.  Wilcox  et  al. 
studied  a  chromium  oxide  producing  premixed  low-pressure  flame  and  observed  an  increase  of 
-180K  1  cm  above  the  burner  upon  addition  of  chromium  precursor.  They  attribute  this  increase  to  a 
shift  in  the  radical  concentrations  towards  equilibrium  by  heterogeneous  radical  recombination 
catalysis  on  chromium  particles.  Lehtinen  and  Zachariah53  predict  that  under  certain  conditions,  par¬ 
ticle  coalescence  will  result  in  a  heat  release  to  the  surrounding  gas. 

The  temperature  increase  was  the  largest  for  the  stoichiometric  undiluted  flame  and  the  smallest  for 
the  rich  flame.  The  temperature  increase  for  the  diluted  flame  was  a  little  less  than  half  that  of  the 
undiluted  flame,  which  is  consistent  with  the  absolute  SiCl4  concentration  being  reduced  by  a  factor 
of  ~0.4  in  the  diluted  flame.  For  the  rich  flame,  the  absolute  temperatures  are  lower,  and  the  tem¬ 
perature  difference  upon  SiCl4  addition  is  relatively  small.  The  lower  temperature  reduces  the  rate  of 
SiCl4  reactions,  which  delays  particle  formation.  At  3  cm  from  the  burner,  4%  of  the  initial  SiCl4 
remains  in  the  rich  flame,  as  opposed  to  <1%  for  the  stoichiometric  flames. 


5.3  OH  concentration 

Addition  of  SiCl4  precursor  depleted  the  OH  concentration  relative  to  the  flames  without  precursor 
for  all  three  basic  flame  configurations.  The  rich  flame  had  by  far  the  most  depletion  on  addition  of 
precursor,  perhaps  because  the  lower  temperatures  prolong  the  processing  of  SiCl4.  The  model  pre- 
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diction  for  the  rich  flame  matches  reasonably  well  in  magnitude  but  not  in  shape.  The  model  predic¬ 
tion  for  the  diluted  flame  is  very  poor,  with  the  model  predicting  much  less  OH  than  is  observed.  We 
did  attempt  to  take  into  account  the  differences  in  efficiencies  of  argon  and  water  as  a  third  body, 
which  did  increase  the  predicted  OH  a  little  but  not  nearly  enough.  The  premixed  burner  reactor  in 
the  CHEMK1N  software  does  not  have  any  provisions  for  including  condensation  or  surface  chemis¬ 
try,  which  could  affect  the  OH  concentration  profile,  especially  in  the  tail  of  the  distribution  where  we 
expect  more  particle  formation  to  be  occurring. 

Zachariah  and  Burgess9  measured  OH  concentration  in  an  atmospheric  pressure  counter-flow  diffu¬ 
sion  flame  doped  with  silane  and  observed  an  order  of  magnitude  depletion  in  the  OH  concentration 
with  respect  to  the  undoped  flame.  They  only  observed  this  depletion  in  the  regions  where  the  pre¬ 
cursor  was  undergoing  pyrolysis  and  particle  formation,  but  not  in  the  particle-laden  zone  after  parti¬ 
cle  formation  was  complete.  They  suggested  that  OH  plays  a  significant  role  in  the  oxidation  of  Si- 
containing  species  and  clusters,  but  doesn’t  react  on  the  highly  hydroxylated  particles.  Hwang  et  al7 
measured  the  OH  concentration  in  an  atmospheric  pressure  co-flow  diffusion  flame  doped  with  SiCl4, 
and  observed  an  order  of  magnitude  reduction  upon  addition  of  precursor  in  the  particle  formation 
zone  but  not  in  other  parts  of  the  flame  (measured  on  the  centerline).  They  attributed  the  OH 
decrease  to  equilibrium  shifts  in  the  OH  concentration  caused  by  consumption  of  O2  and  H20  by 
reaction  with  SiCl4  as  well  as  OH  consumption  by  reaction  with  HC1.  We  observed  suppression  of 
OH  over  our  entire  measurement  range.  We  did  not  observe  the  OH  concentration  profile  with  SiCl4 
approach  its  value  without  SiCl4;  however,  at  low  pressure,  the  particle  formation  zone  is  signifi¬ 
cantly  stretched  compared  to  the  atmospheric  pressure  flames  such  as  those  in  the  experiments 
described  above,  and  all  of  our  measurements  fall  within  the  particle  formation  zone.  The  model, 
which  does  not  take  particle  formation  into  account,  underpredicts  OH  in  the  tail  of  both  the  undiluted 
stoichiometric  and  the  rich  flames,  perhaps  suggesting  that  the  particle  formation  zone  is  beginning  to 
tail  off. 
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6.  Conclusion 


Hydroxyl  concentration  and  temperature  profiles  for  SiCl4  combustion  in  a  low-pressure  premixed 
flat  flame  have  been  obtained  using  LIF.  These  data  are  sensitive  to  the  overall  chemical  kinetics  of 
our  proposed  SiCl4  combustion  mechanism.  In  order  to  obtain  the  best  fit  to  the  observed  data  at  all 
three  flame  configurations,  we  had  to  modify  six  different  rates  from  our  original  estimates.  None  of 
the  modified  rates  are  well  known  in  the  temperature  regime  of  our  flame.  Experimental  study  of 
these  rates  in  the  temperature  range  from  500K  to  1500K  would  allow  us  to  refine  our  SiCl4  combus¬ 
tion  mechanism.  Our  model  did  not  predict  the  correct  OH  profiles  for  all  conditions.  In  particular, 
we  had  the  most  trouble  fitting  the  data  from  the  diluted  flame,  even  though  the  maximum  tempera¬ 
tures  were  similar  to  those  in  the  undiluted  stoichiometric  flame,  indicating  that  our  mechanism  for 
SiCl4  combustion  is  not  complete.  Addition  of  condensation  and  heterogeneous  chemistry  might 
improve  our  understanding  of  this  problem.  While  we  did  not  study  the  effects  of  SiCl4  kinetics  on 
final  particle  characteristics,  we  found  that  even  in  a  premixed  flame,  addition  of  SiCl4  dopant  did 
affect  the  temperature  and  OH  concentration  profile  in  the  flame.  These  profiles  are  determined  by 
the  precursor  kinetics.  Particle  formation  is  very  sensitive  to  temperature  history,  so  precursor  kinet¬ 
ics  could  affect  the  final  particle  characteristics.  The  OH  concentration  was  also  found  to  be  sensitive 
to  the  flame  stoichiometry  and  diluent. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security  programs,  specializing  in 
advanced  military  space  systems.  The  Corporation's  Laboratory  Operations  supports  the  effective  and  timely 
development  and  operation  of  national  security  systems  through  scientific  research  and  the  application  of 
advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical  staffs  wide-ranging  expertise  and 
its  ability  to  stay  abreast  of  new  technological  developments  and  program  support  issues  associated  with  rapidly 
evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  organizations: 

Electronics  and  Photonics  Laboratory:  Microelectronics,  VLSI  reliability,  failure  analysis, 
solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and  CCD 
detector  devices,  data  storage  and  display  technologies;  lasers  and  electro-optics,  solid-state 
laser  design,  micro-optics,  optical  communications,  and  fiber-optic  sensors;  atomic  frequency 
standards,  applied  laser  spectroscopy,  laser  chemistry,  atmospheric  propagation  and  beam 
control,  LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and  evaluation,  battery 
electrochemistry,  battery  testing  and  evaluation. 

Space  Materials  Laboratory:  Evaluation  and  characterizations  of  new  materials  and 
processing  techniques:  metals,  alloys,  ceramics,  polymers,  thin  films,  and  composites; 
development  of  advanced  deposition  processes;  nondestructive  evaluation,  component  failure 
analysis  and  reliability;  structural  mechanics,  fracture  mechanics,  and  stress  corrosion;  analysis 
and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle  fluid 
mechanics,  heat  transfer  and  flight  dynamics;  aerothermodynamics;  chemical  and  electric 
propulsion;  environmental  chemistry;  combustion  processes;  space  environment  effects  on 
materials,  hardening  and  vulnerability  assessment;  contamination,  thermal  and  structural 
control;  lubrication  and  surface  phenomena.  Microelectromechanical  systems  (MEMS)  for 
space  applications;  laser  micromachining;  laser-surface  physical  and  chemical  interactions; 
micropropulsion;  micro-  and  nanosatellite  mission  analysis;  intelligent  microinstruments  for 
monitoring  space  and  launch  system  environments. 

Space  Science  Applications  Laboratory:  Magnetospheric,  auroral  and  cosmic-ray  physics, 
wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric  physics, 
density  and  composition  of  the  upper  atmosphere,  remote  sensing  using  atmospheric  radiation; 
solar  physics,  infrared  astronomy,  infrared  signature  analysis;  infrared  surveillance,  imaging  and 
remote  sensing;  multispectral  and  hyperspectral  sensor  development;  data  analysis  and 
algorithm  development;  applications  of  multispectral  and  hyperspectral  imagery  to  defense,  civil 
space,  commercial,  and  environmental  missions;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  Earth’s  atmosphere,  ionosphere  and  magnetosphere;  effects  of 
electromagnetic  and  particulate  radiations  on  space  systems;  space  instrumentation,  design, 
fabrication  and  test;  environmental  chemistry,  trace  detection;  atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions,  and  radiative  signatures  of 
missile  plumes. 


